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ABSTRACT 

X-ray images of galaxy clusters often display underdense bubbles which are apparently in- 
flated by AGN outflow. I consider the evolution of the magnetic field inside such a bubble, 
using a mixture of analytic and numerical methods. It is found that the field relaxes into an 
equilibrium filling the entire volume of the bubble. The timescale on which this happens de- 
pends critically on the magnetisation and helicity of the outflow as well as on properties of 
the surrounding ICM. If the outflow is strongly magnetised, the magnetic field undergoes re- 
connection on a short timescale, magnetic energy being converted into heat whilst the charac- 
teristic length scale of the field rises; this process stops when a global equilibrium is reached. 
The strength of the equilibrium field is determined by the magnetic helicity injected into the 
bubble by the AGN: if the outflow has a consistent net flux and consequently a large helicity 
then a global equilibrium will be reached on a short timescale, whereas a low-helicity outflow 
results in no global equilibrium being reached and at the time of observation reconnection 
will be ongoing. However, localised flux-tube equilibria will form. If, on the other hand, the 
outflow is very weakly magnetised, no reconnection occurs and the magnetic field inside the 
bubble remains small-scale and passive. These results have implications for the internal com- 
position of the bubbles, their interaction with ICM - in particular to explain how bubbles could 
move a large distance through the ICM without breaking up - as well as for the cooling flow 
problem in general. In addition, reconnection sites in a bubble could be a convenient source 
of energetic particles, circumventing the problem of synchrotron emitters having a shorter 
lifetime than the age of the bubble they inhabit. 

Key words: MHD — galaxies: clusters: general — intergalactic medium — galaxies: jets — 
galaxies: magnetic fields — X-rays: galaxies: clusters 



1 INTRODUCTION 

The gravitational potential wells of galaxy clusters are filled with 
hot (10 7-8 K), hydrostatically - settled gas wh ich emits X-rays via 
thermal bremsstrahlung (e.g. iMolendil 120041) . Many galaxy clus- 
ters, viewed in X rays, display dark caviti es of size ~ 10 kpc at 



various dist a nces fr om the cluster centre jBohringer et al. 1993; 
I Carilli et all 1 19941: iDunn & Fabianl |2004 iMcNamara & NulserJ 
120071 ; iBirzan et al.ll2008h . They are dark because they have a lower 
density than the surrounding intra-cluster medium (ICM), but pre- 
cisely how much less dense is uncertain, except that they are at least 
a factor of three or so less dense. Observationally constraining the 
density is difficult because the line of sight contains also surround- 
ing material; this problem becomes more severe where the bubble 
is at a larger distance from the cluster centre (see lEnBlin & Heinzl 
2002 for details). The bubbles are apparently inflated by an Active 
Galactic Nucleus (AGN) at the cluster centre and then rise buoy- 
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antly through the ICM. In addition, we infer the presence of an 
internal magnetic field and cosmic rays from observed radio syn- 
chrotron emission. 

There is a growing consensus that negative feedback from 
AGN could solve the cooling flow problem: the accretion of gas 
onto a supermassive black hole in the central galaxy releases en- 
ergy to heat the ICM, preventing it from cooling and collapsing 
towards the centre of the cluster. Observationally, there is a strong 
correlation between those clusters which require heating (i.e. have 
a short cooling timescale) and the presence of optical-line em ission 
and radio emission from AGN as well as star f ormation teurnsl 
Il990l : iRaffertv et al.ll2008l : ICavagnolo et al.ll2008l) . The means by 
which this energy might be transferred to the ICM is not yet un- 
derstood, but an interaction between the AGN outflow and the 



surroundings does se e m very likely (e.g. [B riiggen & Kaisei 



2005 



Reynolds et al. 120021: IChurazov et al 1120051 : Briigg en et al 
Brighenti & Mathews 200^). For this reason, it is important to gain 



2002 



some understanding of how AGN-inflated bubbles interact with 
their surroundings. 
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Rising bubbles in a fluid tend to lose their spherical shape 
after rising a distance comparable to their radius. First, bubbles 
tend to flatten while the surrounding medium is flowing around 
them because the material flowing past their sides is moving with 
greater velocity than either in front of or behind them; we know 
from Bernoulli's principle that the pressure at the sides must there- 
fore be lower and so the bubble expands laterally. Then, the bub- 
ble is shredded into many smaller bubbles and eventually becomes 
completely mixed into the surrounding medium. In general there 
is more than one instability responsible for this shredding: the 
Rayleigh-Taylor (R-T) instability appears at the leading edge of the 
bubble where a dense fluid (the ICM) lies above a less dense fluid 
(the bubble) and the Kelvin-Helmholtz (K-H) instability appears 
at the sides of the bubble where there is a discontinuity in veloc- 
ity and density. In the absence of magnetic fields, the growth time 
of the longest wavelength mode (i.e. the bubble radius) of the R-T 
instability is comparable to the time the bubble takes to rise a dis- 
tance equal to its own size. The growth time of the K-H instability is 
likely to be somewhat longer if there is a large density contrast be- 
tween the bubble and its surroundings. However, in many clusters 
we see large bubbles which have risen distances many times greater 
than the ir own size - some mechanism must be inhibiting the insta- 
bilities dJones & De Youn g 2005; R uszkowski et al.ll2007l) . An ob- 
vious candidate is a magnetic field, coherent on the length scale of 
the bubble, either in the ambient medium ('magnetic draping', see 
lLvutikovl[200^ : bursi & Pfrommedl2008h or inside the bubble, or 
both. Alternatively, by analogy with smoke rings it seems plausible 
that there is some purely hydrodynamica l pr ocess responsible - for 
instance Is cannapieco & Bruggenl d2009h and lBriiggen et al.l 12009) 
find that adding a subgrid-turbulence model to hydro simulations 
could encourage the bubble to stay in one piece. Here, I concentrate 
on the magnetic field inside the bubble, and show how an arbitrary 
'turbulent' magnetic field in a new-born bubble could reconnect 
into a large-scale equilibrium and thus provide the necessary rigid- 
ity. This process is similar to that taking place in stars which make a 
transition from convective to non-convective, for instance in proto- 
neutron stars feraithwaite & Spruidl2004lBraithwaitell2008h . 

In sectional look at the process of relaxation to equilibrium, 
finding a relation between the initial magnetic helicity and the equi- 
librium field strength as well as comparing the relevant timescales. 
In section [3] I present numerical simulations of the reconnection 
process, before looking in some detail in section at the structure 
of the equilibria found. In sections [5] and [6] I discuss the results and 
then summarise and conclude. 

2 FORMATION OF EQUILIBRIUM INSIDE A BUBBLE 

In this section I look at the properties of magnetic bubbles and make 
estimates of the relevant timescales and energies. 

Imagine a bubble of hot gas which, after becoming detached 
from the AGN contains an initially turbulent, disordered magnetic 
field. As we shall see below, it is likely that important changes in the 
magnetic field inside the bubble will happen on a short timescale 
compared to the bubble's buoyant rise through the ICM, so we shall 
not consider the interaction with the ICM. Now, generally the field 
left over from the AGN outflow will not be in equilibrium, mean- 
ing that the Lorentz force the magnetic field exerts on the gas is 
not balanced by the gas pressure gradient. Fluid motion results and 
kinetic energy is dissipated by viscosity, which in this low-density 
case is quite high. Eventually the free energy in the magnetic field 



is used up and a stationary equilibrium is reachedQ Before looking 
at the timescale for relaxation to equilibrium in section 12.21 1 first 
calculate the strength of the equilibrium field, making use of the 
conserved quantity of magnetic helicity. 

2.1 Energy and helicity of the equilbrium 

Let us imagine the relaxation to magnetohydrodynamic equilibrium 
of a bubble with radius r, thermal pressure P and density p which 
contains a magnetic field of energy E — VB 2 /Stt, where B is the 
r.m.s. magnetic field inside the bubble and V = 47rr 3 /3 is the vol- 
ume of the bubble. After relaxation, an equilibrium is reached. In 
the following, quantities inside the bubble initially and (finally) at 
equilibrium are marked with the subscripts i and f respectively, and 
quantities outside the bubble - which are assumed not to change 
during the relaxation to equilibrium - have the subscript o. We can 
say the following about the equilibrium state. 

The reconnection destroys magnetic energy on small length 
scales but has little effect on the magnetic helicity, a global quantity 
which is defined as the volume integral of the scalar product of the 
magnetic field with its vector potential H = (l/87r) f A-B dV. 
It can be s hown that in th e case of infinite conductivity, helicity is 
conserved dWoltieJ 19581) . Helicity has units of energy times length 
and so is present more in the larger structures than is the energy - 
and it is approximately conserved during reconnection taking place 
on small scales, a property which has been very useful in many 
conte xts from the laboratory ( Chui & Moffat 19951: iHsu & Bellanl 
l2002h to the solar corona dZhang & Lowll2003h . Therefore: 

H f « (1) 

Consideration of dimensions gives us 

\H f \ =rH,fE f = \ f r f E f , (2) 

where m is the 'helicity length' and Af is a dimensionless param- 
eter of the equilibrium whose value can be determined from sim- 
ulations - in a large-scale equilibrium it should be of order unity. 
Having chosen to put 1/8tt into the definition of helicity, we have 
r H ~ A/B ~ H/E. 

The thermodynamic relations e = P/ (7 — 1) where e is the 
internal energy per unit volume and dU — T dS — P dV give us 

+E f - -Ei = -P (Vf - Vi) (3) 

7 — 1 7 — 1 

where 7 is the ratio of specific heats of the bubble gas and V — 
47rr 3 /3 is the volume of the bubble. It assumed that heat transfer 
between the bubble and its surroundings can be neglected on the 
reconnection timescale, that there is no internal source of heat en- 
ergy, that there is no radiative cooling, so that dS = 0, and that 
the reconnection happens sufficiently slowly that there is always 
pressure balance between the bubble and its surroundings. In fact, 
([3]) is simply an expression of the conservation of enthalpy. Now, 
since the bubble is in pressure equilibrium with the surroundings 
both before and after reconnection, we have 

1 Strictly speaking, the equilibrium should continue to evolve due to the 
finite conductivity of the plasma and the resulting energy loss. We can think 
of the evolving equilibrium as an electric circuit consisting of an induc- 
tance and a resistance, but assuming the standard Spitzer conductivity the 
timescale for decay of the current and magnetic field is ~ 10 33 yrs. 
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where the coefficients 1/3 come from the fact that the magnetic 
field exerts an 'isotropic pressure' P ma g = (1/3)P 2 /Sir, which is 
equivalent to assigning an adiabatic index of 4/3 to the magnetic 
field and using the thermodynamic relation P = (7 — l)e. 

Looking at the equations above, we have four unknowns re- 
garding the final state (rf , Pf , Ef and iff) and four equations 03 to 
(|4]). [The first equality of (|4} doesn't count, as it relates only quan- 
tities in the initial conditions.] So it is now possible to solve for the 
final state once we know the value of Af . 

Now, note that if the plasma is relativistic we have 7 = 4/3 
and Vf = Vi since both plasma and magnetic field have the same 
ratio of pressure to energy density, so converting energy from one 
form to the other has no effect on the total pressure. If on the other 
hand the fluid is non-relativistic and monatomic, i.e. 7 = 5/3, 
we expect a fractional increase in the bubble's volume; the great- 
est increase in volume possible occurs if the bubble contains only 
magnetic energy and no thermal energy at the beginning and where 
all of this energy is converted to thermal, i.e. where Pi — and 
Hi — 0, and it is easily found that the fractional increase in volume 
Vf/Vi — 8/5. Making the approximation that Vf w Vu we can 
retrieve from the above set of equations the following relation 



B[ 
B? 



A f ' 



(5) 



where B is the r.m.s. field strength in the bubble. Note that the 
final magnetic energy and volume of the bubble are not dependent 
on the gas density; the latter will certainly affect the time taken to 
reach equilibrium, and may or may not affect which equilibrium is 
reached from a given initial magnetic field, but once an equilibrium 
is reached the density no longer has any effect. I return to the effect 
of density in section UA\ If we need to know the density in the 
equilibrium bubble we have pf/pi = Vi/Vf , but in any case we see 
from above that the density can drop by no more than a factor 5/8. 



2.2 Timescales 

Imagine an initially stationary fluid of uniform gas pressure con- 
taining an arbitrary magnetic field. Looking at the momentum 
equation and comparing the sizes of the terms, we see that fluid 
velocities comparable to the Alfven speed will be induced, and so 
the relaxation to equilibrium must take place at this speed. Gen- 
erally the equilibrium will be topologically different from the ini- 
tial conditions and so reconnection of field lines will be required; 
this reconnection should also take place at the Alfven speed, re- 
gardless of the mechanism invoked and its microphysics. Studies 
of reconnection in various contexts confirm that reconnection does 
take place at roughly that speed, or rather, somewhat less, say avA 
where a has a value of around 0.1 dElsner & Lambll 19841 ; llkhsanovl 
1200 ll) . This means that the time required for reconnection to occur 
across a structure of size / rec is T rec ~ Zrec/ ((xva). 

Now, imagine that a bubble of radius r contains a magnetic 
field structured on some length scale /. Initially, / = h and re- 
connection proceeds on this length scale but as the field relaxes 
towards equilibrium, the length scale of the magnetic field grows 
until eventually / « r. This means that the magnetic field should 
initially evolve on a timescale of T rec ~ 1\/(olva) but that as it ap- 
proaches equilibrium T rec ~ r /(avx)- At the same time, magnetic 
energy has been dissipated and va has fallen, by a large factor if 
k <C r. In fact, the magnetic energy should evolve according to 



dln_E 

dt 



1 

'Tree 



avA 



(6) 



Clearly the magnetic field evolves slowest when it is closest to a 
global equilibrium, therefore the time taken to relax to the equi- 
librium is approximately equal to the reconnection timescale with 



"7~relax 



avA 



7.1x10" 



B 



aB 



yr, 



6E 



lOkpc 



10 5 m p cm 



(7) 



(8) 



\20pG 

where va, B and E are the Alfven speed, magnetic field and energy 
at equilibrium. 

These bubbles are embedded in the hydrostatically- settled in- 
tracluster medium and, being less dense, they rise upwards through 
it. We can calculate the buoyant rise velocity and the associated 
timescale, which we define as the time taken for the bubble to move 
a distance equal to its radius, as this is the timescale on which the 
bubble might be disrupted by instabilities at its surface. Equating 
the buoyant force to the drag force, we have 



gV(po - p) = -p u SC d 



(9) 



where g is gravity, p Q is the density of the external medium, u is the 
terminal velocity, S is the cross-sectional area of the bubble and Cd 
is the drag coefficient, which has a value of aroun d 0.5 for a solid 
sphere at Reynolds numbers of ~ 10 5 . However. IChurazov et al.l 
(1200 II) find a higher effective value in the context under considera- 
tion here, owing to extra energy loss from the excitation of internal 
gravity waves in the ambient medium; below, we shall adopt their 
value of Cd = 0.75. Also, we can express gravitational accelera- 
tion g in terms of the Keplerian velocity ^Ke P of a circular orbit at 
distance R from the centre of the cluster, g = v^ ep /R. Assuming 
the bubble is spherical, the terminal velocity u is given by 



8 p - p r 



;^Kep- 



(10) 



3Ca po R 

Note that the rise velocity is comparable to the Keplerian veloc- 
ity, which in turn must be comparable to the sound speed in the 
hydrostatically- supported surrounding medium c s , ; the relation 
between the two is v^ ep /R = c^ /(j Hp) where Hp and 70 
are the pressure scale height and adiabatic index (= 5/3) of the 
surrounding gas. However, the motion of the bubble is not likely to 
be fast enough that significant energy is dissipated in shocks. The 
rise timescale is 



Trise — 



SCd 



R 



VKep 



5xl0 6 f- 



^Ke P 



lOkpcy V 1000 km s 



=0- 



(11) 



yr, (12) 



taking p /(p — p) ~ 1, which seems justified by the X-ray ob- 
servations. As a check, we can also calculate the time taken for the 
bubble to reach this terminal velocity u by looking at the accelera- 
tion a from rest: 



7~acc — — 



up 



8 r 
3Cd R 



R 



(13) 



a d(Po - p) \ '6Ud ti J VKep 

where p is some density which accounts for the inertia of the bub- 
ble and of the surrounding gas; it is comparable to p Q . Compar- 
ing this to the rise timescale at terminal velocity r r i se we find that 
Tacc/T r i S e ~ 8/3Cd which is roughly equal to unity. We therefore 
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make the approximation that the bubble always moves at its termi- 
nal velocity. 

We can now compare the relaxation timescale r re iax to the rise 
timescale. Dividing by dTTb we find that 




On the second line above, the ratio has been expressed as the prod- 
uct of two factors of order unity in addition to a pressure ratio, a 
density ratio, and a factor of 1/a. The quantity P ma g = B 2 /247T 
is the isotropic magnetic pressure (a quantity which appears again 
in section|4]) inside the bubble. It is these ratios P /Pma g and p/p G 
which are most uncertain; it seems likely that the former is rather 
high and that the latter is rather low. In the third line some likely 
numerical values are given which produce a ratio of timescales of 
about 1. 

In light of these estimates, it is natural to ask what the mag- 
netic field should look like if we observe it before it has reached a 
global equilibrium. Observationally, there are the following possi- 
bilities: 

(i) The initial field is weak and no significant reconnection oc- 
curs. The dominant length scale in the field / remains at its initial 
value li where T rec ~ 1\/(olva) > i~ age where r ag e is the age of the 
bubble. The field evolves passively in response the the interaction 
between the bubble and its surroundings. 

(ii) The field is stronger and reconnection proceeds initially on 
a short timescale, but because of low helicity the equilibrium field 
strength is low and the Alfven speed drops by a large factor as 
reconnection proceeds. At the time of observation, the field is still 
structured on scales small compared to the size of the bubble: li < 
I < r. If it were possible to measure / and va, we would find that 

l/(0LVA) ~ Tage- 

(iii) The helicity of the field is large so that the equilibrium en- 
ergy as calculated in © is also large. Reconnection proceeds and 
the length scale of the field grows without the Alfven speed becom- 
ing very low. A global equilibrium is reached (/ w r) and we would 
measure that r j (old a) < r age - 

Given the large uncertainly in the parameters, it is not at all 
clear whether the field inside the bubble should have time to reor- 
ganise itself into an equilibrium before the bubble rises and is dis- 
rupted into a 'mushroom-cloud' shape. The magnetic field strength 
in observed bubbles could well be somewhat greater than 20pG, or 
the density could be much lower - indeed the material in the bubble 
could be a pair plasma instead of ionised hydrogen - both leading 
to effective relaxation on a short timescale and possible stabilisa- 
tion of a bubble against shredding instabilities, but it is impossible 
at present to say for certain whether the ratio of the two timescales 
Q3} in observed systems is less than, equal to, or greater than unity, 
and it is possible that all three regimes exist in different bubble sys- 
tems, given the diversity in observational properties between dif- 
ferent galaxy clusters. 



2.3 Input from AGN outflow 

We have seen above how the evolution of a bubble depends on the 
magnetic helicity of the field it contains, as well as the field strength 
- it would be useful therefore to look at how much magnetic energy 
and helicity an AGN is likely to contribute. 

The bubbles are presumably inflated by some jet or other out- 
flow from a system consisting of a supermassive black hole and 
an accretion disc. I assume here that the bubble is inflated by a 
magnetocentrifugally-accelerated j et from an accretion disc with 
a net- flux magnetic field (see e.g. Bland ford & Pavndll98ilMoiil 
2009, and refs therein). The net poloidal flux comes from the ac- 
creted material and may therefore change in time, but assuming that 
at any one instant the disc is threaded by magnetic field of a par- 
ticular direction, the material in the jet will also contain a poloidal 
field in this direction. The direction of the toroidal field in the jet 
will also depend on the direction of rotation of the disc, such that 
the toroidal field will be in the same sense as the rotation of the disc 
if the poloidal field is directed into the disc and it will be in the op- 
posite sense to the rotation if the poloidal field is directed out of the 
disc. Since we can think of helicity as the product of poloidal and 
toroidal fluxes, we see that the two jets, and therefore bubbles, on 
either side of the disc will have equal and opposite magnetic helic- 
ities. This is obviously convenient because the disc does not need 
to produce any helicity itself; rather, it merely transports helicity 
between the two hemispheres. 

The following is a rudimentary estimation of the magnetic 
energy and helicity in a bubble. Although AGN jets are thought 
to have bulk Lorentz factors ~ 10, the calculation here does not 
take relativistic effects into account: while there is no certainty that 
the bubble is in fact inflated by a relativistic outflow - it is likely 
that the majority of mass and magnetic flux comes from a non- 
relativistic disc wind - there seems little purpose in conducting a 
relativistic generalisation at this stage. 

The system can be characterised by the half- width of the out- 
flow at the Alfven surface r- } as well as the gas pressure Pj, the den- 
sity pj, the flow speed Vj and the magnetic field Pj at the Alfven 
surface. At the Alfven surface, the poloidal and toroidal compo- 
nents of the field are roughly equal, and noting that helicity is the 
product of poloidal and toroidal fluxes we may write down the fol- 
lowing expression for the helicity crossing the Alfven surface per 
unit time: 

H t ~ $ P oi$tor,t - Bfrfvi (15) 

where the subscript t denotes a time derivative. All factors of or- 
der unity are dropped. The helicity per unit mass injected into the 
bubble is this quantity divided by the mass injection per unit time 
Pjrfvj, so that the total helicity injected into a bubble of radius rh 
and density pb (using the findings of section IXTl that n « n and 
Pi w pf and using the subscript b for 'bubble') is 

H~Bf ri rld. (16) 
Pi 

If there is no net flux through the disc, or if the direction of the 
net flux changes polarity from time to time, the bubble will con- 
tain very much less helicity than this estimate, leading to an even 
weaker field inside the bubble once reconnection has taken place. 
Incidentally, the jet or outflow will presumably go through a shock 
as it passes into the bubble - I make the approximation here that 
the shock does not dissipate significant magnetic energy or helic- 
ity. Imagining that the bubble inflates without any conversion of 
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magnetic energy to heat or vice versa, so that the magnetic field 
expands adiabatically0 



4/3 



(17) 



where B{ is the field strength in the bubble after inflation but before 
reconnection has begun. Of course, in reality reconnection will be- 
gin whilst the bubble is still inflating, so that the prediction here of 
whether the field strength is high enough for reconnection to pro- 
ceed at all can be considered rather conservative. 

Now, once inflated the bubble field relaxes into an equilibrium 
with a dimensionless helicity length Af w 1 whose field strength 
Bf is given by 



Hr7 



rb Pi ' 



(18) 



where d~i~6T) has been used. 

The bubble moves away from its 'mooring' when the expan- 
sion speed drb / dt slows down to approximately the buoyant rise 
speed u, which using dTTTb and considering the mass flux into the 
bubble pjVjrf, gives 



Pb vi 

Pi VKep \Vh 



(19) 



which can be used to eliminate the ratio pb / Pj in (TTvT) and (QjO), re- 
placing it with the ratio Vj / VKep whose value is better constrained. 
Dividing (fl"8l) by ( fTTt and using $5$ and dT9t , taking Af « 1, gives 



x^ri el 

r h \pb 



1/3 



rb Vi J 



1/3 



(20) 



In a typical AGN/bubble system, we might expect r- } « 10 cm, 
rb ~ 10 22 cm, V) ~ c and VKep ~ c/300, which gives Ai ~ 10 -3 . 
Therefore we expect the energy of the magnetic field to fall by a 
large fraction during relaxation to equilibrium. The length scale of 
the field in the bubble before reconnection begins is given by 



1/3 



/ \ 1/3 / \ 2/3 



(21) 



from consideration of the expansion of material from the outflow 
and using ( TT9t . This is consistent with H ~ Erfi which is a special 
case of the general inequality H < El (or alternatively A < l/r) in 
the case where the twist is consistently in one direction, like a box 
full of right-handed screws. 

Defining a magnetisation parameter 



Pic 2 



(22) 



it is possible to express the timescale ratio (TUT) in the following 
way (again, ignoring factors of order unity): 

1/3 

(23) 



"7~rela 



1 /^Kep\ 7 / 6 / Tb\ 

sfm V C J V r j / 



2 Strictly speaking, I have also assumed that the expansion is isotropic. 
Whilst this is clearly not the case in a jet, I am assuming that the expansion 
on the other side of the shock cancels this effect. The justification for this 
is the tendency of any magnetised volume to adjust its shape until the three 
components of the field are roughly equal since this represents the minimum 
energy (see section ITU . Since the three components are roughly equal at 
the Alfven surface, the total expansion must be the same in each dimension. 



If the outflow is modestly relativistic with m ~ 10, we have 
a timescale ratio of around 10 with these parameters. Given the 
approximate nature of this calculation and the uncertainty in the 
parameters, we cannot therefore say with confidence whether the 
magnetic field should reorganise into global equilibrium before the 
bubble moves far. However, we can estimate the initial reconnec- 
tion timescale T rec ,i: 



(h/avA,f) 



Ti VKep 



1/2 



(24) 



7"relax (r h /av A ,i) r h H y \r h Vj 

using j2Tb . and J20b . Since this has a value < 10 -4 and there- 
fore Tree , i / Tr ise < 10 , we should certainly expect to see re- 
connection in progress if an equilibrium has not already formed. 
In terms of the three cases described in section 12.21 case (i) looks 
very unlikely unless the outflow is extremely weakly magnetised, 
in which case it must be driven by some non-magnetic mechanism. 



3 NUMERICAL SIMULATIONS 

In this section, simulations of the relaxation of a turbulent field into 
equilibrium are described. 

3.1 Numerical scheme 



The code used is the STAGGER CODE dNordlund & Galsgaardl 
I1995L iGudiksen & Nordlundll2005h . a high-order finite-difference 
Cartesian MHD code which uses a 'hyper-diffusion' scheme, a 
system whereby diffusivities are scaled with the length scales 
present so that badly resolved structure near the Nyquist spatial 
frequency is damped whilst preserving well-resolved structure on 
longer length scales. This, and the high-order spatial interpolation 
and derivatives (sixth order) and time- stepping (third order) in- 
crease efficiency by giving a low effective diffusivity at modest res- 
olution (144 3 here). The code includes Ohmic and well as thermal 
and kinetic diffusion, which are kept at a low level. The code uses 
Cartesian coordinates with periodic boundaries, although the com- 
putational box is made sufficiently large that nothing significant is 
happening at the boundaries. 

3.2 Numerical setup and initial conditions 

The bubble is modelled as a sphere of initial radius n containing 
hot gas (with an ideal gas equation of state and 7 = 5/3) and a 
turbulent magnetic field, surrounded by a cooler, unmagnetised gas 
with the same equation of state. 

The size of the computational box is chosen so that there are 
no issues with the bubble material expanding across a boundary 
and back into itself; it is found that using a computational box of 
size 6n suffices. In these calculations, there is no gravity, and the 
ambient medium has a uniform density as we are interested primar- 
ily in processes happening on timescales shorter than any buoyancy 
timescale. 

At this stage, we assume the ambient medium is unmagne- 
tised. While this is strictly speaking probably not the case with 
the radio bubbles observed, it is safe to assume that the Alfven 
speed inside the bubble is significantly greater than that outside 
and that therefore during the timescales of interest here, the ambi- 
ent medium will not evolve. Also, it is possible that the ICM was 
originally not significantly magnetised and that it became magne- 
tised (and chemically enriched) by galactic mass ejection, whether 
from AGN or from stars. 
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The bubble is given an initially random magnetic field which 
contains energy at a range of length scales: the minimum wavenum- 
ber in the initial field IS /Cmin = 27r/n, i.e. the largest length scale 
present is equal to the bubble radius, and the energy declines to 
higher wavenumbers as E(k)5k oc k~ 5 ^ 3 (for want of anything 
better). 

From the considerations in section 12.11 we have the follow- 
ing degrees of freedom in the initial conditions once we have ar- 
bitrarily fixed n and P Q and chosen 7: an initial helicity parame- 
ter Ai = Hi/{r\E\), the density ratio pi/p Q and the pressure ra- 
tio Pmag/i^o- In the simulations described below, a range of these 
parameters is explored. In section 13.31 simulations are run with 
initial ratios pi/p Q — 1/10 (consistent with observations) and 
Pmag/Po = 1/2 and the effect of the magnetic helicity is exam- 
ined. Later, other parameters are looked at. 



3.3 Dependence on initial helicity 

In this section, a set of simulations with different values of the di- 
mensionless helicity parameter Ai is presented. As the field evolves 
on the dynamical Alfven timescale, we see how the magnetic en- 
ergy E, helicity H and other parameters change. The evolution of 
energy and helicity in these simulations is plotted in figs. [TJ [2 and 
[3] The simulations are labelled with letters a-g in the figures. All 
simulations were run for the same number of timesteps, which cor- 
responds approximately to the same number (about 600) of sound- 
crossing times t\I c where c Q is the sound speed in the external 
medium. Clearly, whilst the energy falls dramatically at the begin- 
ning of each run, the helicity falls much more modestly. In fact at 
first helicity does not seem to fall at all while the energy drops by 
a factor of ten or more. Then, it is found that in some cases, a sim- 
ple equilibrium is reached after a short time whereas in other cases 
the field loses more energy and any equilibrium is often more com- 
plicated in shape - helicity is the determining factor, as expected. 
In cases with high helicity, less energy is lost and an equilibrium 
is reached more quickly. It seems though that in all cases some 
equilibrium is eventually reached, after a number of Alfven cross- 
ing times (see fig. O, the difference being that in some runs the 
Alfven crossing timescale becomes very long before equilibrium is 
reached. Also, we see that at equilibrium, all runs lie close to the 
straight dashed line in fig.[3j \H\ = OAnE, which reflects the fact 
that all equilibria have comparable helicity lengths |Af \n ~ 0.4n, 
defined in ©. 

In addition, note that in fig.[3]it can be seen that once an equi- 
librium is reached, energy and helicity fall together. This is an effect 
of finite conductivity; energy should still fall somewhat faster than 
helicity since the length scale of the equilibrium rises; the equilib- 
rium 'spreads out'. To be more precise, we see from integrating 
over volume the zero-velocity diffusion equation d~B/dt = 77V 2 B 
and using Gauss' theorem to equate the right-hand side to zero that 
BV ~ Br 3 ~ const, so that E oc r -3 . Now, since H ~ rE we 
have H oc E 2/s as the equilibrium diffuses (illustrated in fig.O. 
Also note that flux <£ oc if 1/2 . 

As the magnetic field evolves, it was found in section [27TI that 
the volume of the bubble increases, and in this case it should in- 
crease by a factor 13/10 if the magnetic energy is largely converted 
into thermal. In the simulations, some increase in bubble volume 
is seen, although it is not possible to draw firm conclusions as to 
whether the volume increase is just that predicted or whether there 
is an additional increase from mixing of the ambient medium into 
the bubble. This topic will be explored in a forthcoming publica- 
tion. 




500 



t / 



Figure 1. Log magnetic energy against time, the former in units of P rf 
and the latter in units of the Alfven timescale at t = 0, i.e. tao = 
ri\/piVi/(2Ei), for a set of simulations with various different values of 
the helicity parameter Aj. Where the lines are dotted the field is non- 
equilibrium, solid lines signify some kind of non-simple equilibrium and 
thick solid lines signify that the field is in a simple axisymmetric equilib- 
rium. The initial energy is the same in all cases (E\ = 2tt in these units) 
but clearly the eventual energy varies between the simulations. 




Figure 2. As fig. [T] but the time axis has been normalised to give the true 
number of Alfven timescales which has elapsed. The Alfven timescale in- 
creases as the magnetic energy falls so that the time unit in the plot is 
increasing to the right. In this plot we can see that equilibrium is indeed 
reached on the order of ~ 10 ta- 



In fig. [4] the evolution of the magnetic field in one of these 
simulations (the one marked 'a') is illustrated. The field has a rela- 
tively high helicity (Ai = 0.032) and the field quickly evolves into 
a simple torus equilibrium, via a figure-of-eight shaped configura- 
tion consisting of a twisted flux tube twisted around itself. 

In fig. [5] the evolution of the magnetic field in another sim- 
ulation (f) is illustrated, this time with lower helicity parameter 
Ai = 0.0013. An intermediate state is reached which consists of 
two torus- shaped fields connected by two flux tubes; however the 
two tori are pulled together again by the tension in the tubes which 
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Figure 3. Log magnetic energy against log helicity for the same set of sim- 
ulations as in fig. [T] Energy and helicity are plotted in units of P rf and 
P Q rf. As in fig.H where the lines are dotted the field is non-equilibrium, 
solid lines signify some kind of non-simple equilibrium and thick solid lines 
signify that the field is in a simple axisymmetric equilibrium. The stars show 
the position at the end of the run. Helicity falls only a little during the ap- 
proach to equilibrium, and there is a clear correlation between helicity and 
the energy at which the equilibrium is reached - the straight dashed line 
illustrates the relation H = QAr\E. The solid lines in the lower-right show 
the gradients at which the field should move during reconnection (which 
conserves helicity) and during pure Ohmic dissipation (H oc E 2 / 3 ). 

join them, and a current sheet forms at the interface. Eventually the 
two tori become one, significantly weaker torus. 

It appears that the most basic equilibrium is a 'ball of string' 
twisted torus shape. All equilibria consist of twisted flux tubes ar- 
ranged in some pattern, the simple torus being a special case where 
the flux tube makes a circle. Tubes may be twisted in either sense, 
corresponding to positive and negative magnetic helicity. Magnetic 
fields with greater helicity tend to evolve directly into simpler equi- 
libria; conversely when the helicity is very small the energy drops 
by a large factor, reducing the Alfven speed to such an extent that 
continued simulation of the evolution becomes impossible. It seems 
plausible that all fields eventually evolve into a simple torus con- 
figuration, the important question being whether this happens on a 
sufficiently short timescale. 

3.4 Density of the bubble 

We now explore the effect of the bubble density. To do this, simula- 
tions were run with the following values of the density ratio p\ / p Q 
parameter: 0.1 (as above), 1 and 10. Of course, the former (0.1) 
is the only ratio consistent with the observations but it is informa- 
tive to look at other values. In fig.[6]three simulations are compared 
which have identical initial conditions except for p\/ p Q \ their be- 
haviour is quite different. 

The main difference between the simulations is that the denser 
bubbles move into the surroundings more easily, become more non- 
spherical and are more prone to breaking up, even though the field 
strength is the same and the Alfven speed is lower. In the fig- 
ure we see that in the low-density case (left column) the bubble 
becomes somewhat distorted but then returns to a more spherical 
shape, forming a simple torus. In the run with pi/p Q = 1 (middle 
column) different parts of the bubble move away from each other 



and several small torus shapes are formed, connected by weak flux 
tubes. 

To illustrate how a bubble can break up, field lines of a run 
with pi/po = 1 are plotted in fig. [7] This run has very low he- 
licity Ai = 0.00025, and the magnetic field is identical to run 'g' 
described in the previous section. The bubbles with higher helicity 
always form a simple torus-shaped equilibrium regardless of the 
density. 

This dependence on the bubble density can be understood in 
the following way. During reconnection, material inside the bubble 
is moving around with velocity comparable to the Alfven speed, 
so that the bubble will not stay spherical for very long. The kinetic 
energy density of the plasma is comparable to the magnetic energy 
density B 2 /8tt and its momentum per unit volume is B\J p/An. A 
denser (and colder) bubble has more momentum and can penetrate 
the surrounding medium more easily. Alternatively, one can think 
of the distance over which a projectile slows down via aerodynamic 
drag: it is comparable to the distance over which it has to push its 
own mass out of the way, which is obviously further if it is more 
dense. Since in reality we know that bubbles have a low density, 
i.e. p\/ p < 1/3 and probably even lower, we should expect only a 
modest deformation of the bubble from the Alfvenic motions inside 
them. 



3.5 Effect of resolution 

It is often useful when employing numerical methods to examine 
the effect of resolution. To this end, simulations were run at double 
the spatial resolution (i.e. half the grid-spacing Ax) used in the 
simulations described in the previous sections. To avoid excessive 
computational cost, the computational box was made smaller: the 
domain had sides L = 4.5n, as opposed to 6n as used previously; 
in the particular simulations run, it was checked that the boundaries 
did not cause significant problems. To be sure of separating the 
effects, some low resolution simulations were run with L — 4.5ri, 
so that grid spacings of Ax = n /24 (as used previously) and n / 48 
could be compared. In fig.[8]the evolution of the energy and helicity 
of the magnetic field in two otherwise identical sets of simulations 
are presented. 

It can be seen in the figure that helicity conservation improves 
at higher resolution. Physically, higher resolution is equivalent to 
lower diffusivity (kinetic, thermal and magnetic), and we would 
expect that this causes the thickness of reconnection sheets to fall. 
This is because the reconnection speed is set by the Alfven speed 
and whatever microscopic (diffusive) processes are involved there- 
fore have to work at a given rate; lower diffusivity means that the 
process has to take place on a smaller length scal^E This in turn 
means that less helicity is destroyed in the reconnection zones, 
since helicity has units of energy x length. This validates the as- 
sumption of approximate helicity conservation 0. Finally, note 
that increased resolution causes the equilibrium, once formed, to 
evolve diffusively (downwards and to the left in the fig. [8} more 
slowly. 



3 That the relaxation of a magnetic field into equilibrium involves current 
sheets at all can be seen from the same argument. A rigorous proof of the 
formation of current sheets in three-dimensional magnetic relaxation does 
not exist, however: iGru zinov ('2009j) discusses a proof of current sheet for- 
mation in two dimensions. 
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Figure 4. Time sequence from the highest helicity run (a) of those in figs.[T]to[3] Left: initial conditions. Centre: figure-of-eight configuration. The field can 
be thought of as consisting of a single flux tube twisted around itself in the opposite direction to its internal twist. Right: later, once a (roughly) axisymmetric 
equilibrium is reached: the tube has become a circle. Such an equilibrium can have either negative or positive helicity, corresponding to a clockwise or 
anticlockwise twist in the field lines. 



Figure 5. Time sequence from the second- 
lowest-helicity run (f) , starting from top- 
middle, then top-right, then bottom left, etc. 
Two tori form, connected by two thin flux tubes. 
Bottom-left: the flat surface is coloured accord- 
ing to field strength: red is weak, green/blue is 
strong. Bottom-centre: the tension in the flux 
tubes pulls the tori together, and a small region 
of high current density is visible, represented by 
green volume rendering. Bottom-right: the two 
tori have merged. They had opposite sign helic- 
ity, reflected in the opposite sense of the twist; 
the resultant torus, which has much lower en- 
ergy than either of the two original tori, has the 
same sense twist as the larger of the two. 




4 STRUCTURE OF THE EQUILIBRIA 

In an MHD equilibrium in a medium without gravitational force, 
the Lorentz force (l/47r)(V x B) x B, which is perpendicular to 
the magnetic field, is balanced only by the pressure gradient force 
—VP, so we can see that the magnetic field lines lie on surfaces of 
constant pressure. Since pressure is a scalar field, the equilibrium 
must consist of nested magnetic surfaces. This is clearly inacces- 



sible from the turbulent initial conditions considered here without 
some topological reorganisation, i.e. reconnection, hence some loss 
of magnetic energy is inevitable. 

We see from the simulations that the basic building block from 
which the various equilibria are made is the twisted flux tube; in 
the simplest case which I describe in more detail below, there is 
a single circular-looped twisted flux tube, and in other cases the 
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Figure 6. Top row: contemporaneous 
snapshots from three otherwise iden- 
tical simulations with pi/p = 0.1, 
1 and 10 (left to right respectively). 
The 0.1 case is run 'e' in the previ- 
ous section. Bottom row: the two bub- 
bles with pi/po = 0.1 and 1 at a later 
time; in the highest density case the 
bubble spills over the boundaries of 
the computational box and no second 
snapshot is plotted. 




twisted flux tubes are arranged in more complex patterns. I shall 
now explore the properties of these tubes. 

4.1 Twisted flux tubes 

The magnetic configurations reached from lower-helicity initial 
conditions in general consist of intertwined, branched, twisted flux 



tubes. The tension in a flux tube is proportional to 2B^ — B\, 
where B ax and are the axial (parallel to tube axis) and hoop 
(perpendicular to tube axis) field components; see the Appendix 
for a proof. This means that a flux tube with only a modest twist, 
i.e. with Bh <C B ax , will have a tension, tending to shorten the 
tube until 2B^ — B\. The tubes connecting the two tori in fig. [5] 
are connected by such untwisted flux tubes, with the result that the 
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Figure 8. Log energy against log helicity for two sets of simulations per- 
formed at standard resolution Ax = rj/24 (dotted lines) and double that 
resolution (solid lines). It is very clear that with improving resolution, he- 
licity is better conserved during relaxation to equilibrium. 

two tori are pulled towards each other and merge. In contrast, the 
flux tubes visible in the top-centre frame of fig. [6] are twisted, with 
the result that the structures are not pulled towards one another and 
at a later time (lower-right frame) are still at similar distances. 

A twisted flux tube will in general be subject to an interchange 
instability of low azimuth al wavenumber JShafranovlll956l : iTavled 
Il95 ilKruskal et al.lll958h . It is found that in the limit of low diffu- 
sivities, the tube becomes unstable to a m — 1 'kink' mode when 
the number of field-line windings around the tube exceeds unity, 
i.e. when 

->2tt^. (25) 
a B h 

In light of the zero-tension condition 2B 2 ^ = B^, this in effect sets 
an upper limit on the length of a flux tube. The instability will result 
in reconnection, loss of magnetic energy and hoop flux <Iv As $h 
falls, the tube will tend to contract and become wider, thus restoring 
equilibrium. The instability will be quenched as the left-hand- side 
of ( |25t falls and becomes equal to the right-hand- side. 

The tendency of a long flux tube to become shorter in this way 
is in some sense a manifestation of the tendency of a flux tube of a 
given helicity to reduce its energy (although it is not clear whether 
the instability just described does in fact conserve helicity). Now, 
if we allow <J> ax and $h to change whilst holding their product H 
constant (as might be the situation during initial reconnection) then 
we have (ignoring some factors of order unity) 

H ~ $ ax $ h « 7ra 3 lB^B h « aVB 2 , (26) 

the last step following from the equilibrium condition (IA121 ). The 
energy of the field E ~ VB 2 , so that 



This means that the lowest energy state of a flux tube with a given 
volume and helicity is / — > 0, a — > oo. The reason that short, 
fat flux tubes do not actually appear in the simulations presumably 
has something to do with the boundaries at the ends of the tubes 
and with limits on the flux <J> ax through these boundaries. One can 
imagine various kinds of end-boundaries for a flux tube, the sim- 
plest being periodic, or in other words, the tube is connected to 



itself in a loop. In this case, it is clear that the length / ^ 27ia, and 
although of course the assumption of a straight tube with circular 
cross-section made above will no longer be accurate, one would 
certainly expect to find a tendency for tube loops to contract as 
much as possible and reach / « 2ita (via kink instability once a 
tube has formed, or otherwise); this case where a flux tube loop is a 
circle is discussed in the next section. Also possible is that the tube 
merges with other tubes (some of which might resemble simple cir- 
cular loops), perhaps in such a way that the ends of the tube may be 
considered anchored and that <J> ax is fixed. This kind of tube can be 
seen in the top-centre and top-right frames of fig. [6] where various 
flux tubes are attached to each other in various ways. 

One flux-tube arrangement which often makes an appearance 
in the simulations is the figure-of-eight and variations thereof. In 
its simplest form it consists of a single flux tube wrapped around 
itself in a double-helix fashion. Its total length is apparently greater 
than the kink-instability limit given by (|25l > and dA12K which could 
be allowed because this configuration is in fact the end- state of the 
instability or because there is some degree of branching, with field 
lines passing between neighbouring sections of the tube. Whilst 
this question and the detailed structure of this configuration is left 
for future investigation, it is possible to understand intuitively why 
the tube should twist around itself in this manner. Imagine starting 
with a straight flux tube with given axial flux but with zero twist, 
i.e. Bh = 0. Now rotate one end of the tube by an angle a, and 
note that Bh oc a and that the hoop-component energy Eh oc a 2 
so that the torque required to produce an additional rotation 5a 
is proportional to a. Now imagine two parallel tubes - both with 
Bh = - touching along their length and joined somehow at the 
ends (so strictly speaking, one flux tube). Now rotating one end of 
the configuration by an angle (3 about the axis created by the line 
of contact of the two tubes will produce a new field component 
perpendicular to that axis and it is easily verified that the torque 
required to produce an additional rotation is again proportional to 
the angle f3. Finally, imagine a configuration where the two tubes 
have an initial Bh as if they had been twisted by an angle a; the 
tubes will exert a torque on the ends in such a way that a helical 
rotation will be produced in the opposite direction from the internal 
twist of the field lines inside the tubes, until the two torques are 
in balance. Since the tubes are touching and therefore the radii of 
rotation are comparable, in the end state a ~ /3. This approximate 
relation is consistent with the simulations. 

It is not clear exactly how the transition is made from the 
double-helix configuration to the simple torus, although it is com- 
mon in the simulations (see fig-SJ). However, there should be some 
drop in hoop flux in relation to axial flux. A deeper investigation of 
this point is left to the future. 



4.2 Simple axisymmetric equilibrium 

The most basic self-contained equilibrium is that shown in the 
right-hand frame of fig-H] This configuration is approximately sym- 
metrical about some axisQ In an axisymmetric equilibrium the az- 
imuthal co mponent of th e Lorentz force must vanish, leading to the 
condition jMestelll 19611) : 

Bp • V(^£ ) = 0, (28) 



4 A way of predicting the direction of this axis from the initial conditions 
has so far escaped investigation. 
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Figure 9. Contour plot of ip in the meridional (vj,z) plane in the equilib- 
rium illustrated on the right-hand- side of fig.H] From the definition of ip, the 
contours are also poloidal field lines, whose direction is indicated by the ar- 
rows. The contours are equally spaced, meaning that there is equal poloidal 
flux between each of the plotted field lines. [Note that one consequence of 
this is that the density of lines plotted vanishes towards the axis, although 
the poloidal field is strong in this central region.] In addition, the left and 
right hand sides contain contour plots of pressure difference P — P and 
of F = tzjB^ respectively, represented by red and blue shading. It is clear 
that P = P(ip) and F = F(ip) as predicted above. 



are the poloidal (meridional) and toroidal (az- 
imuthal) components of the magnetic field and w is the cylindri- 
cal radius. Everything can be described in terms of a flux function 
ij)(vj,z) which is the product of the cylindrical radius and the az- 
imuthal component of the vector potential A, i.e. ip = wA$ where 
Bp = V x (A^e^), so that 



zuB z 



dip 
P = P(ip) 



and 
and 



voB- 



_chp_ 
dz' J 

wB <l> = F{^). 



(29) 
(30) 



The condition that the Lorentz and pressure gradient forces ba lance 
can now be expressed (Grad & Rubin 1958; Shafranov 1966) as: 



A * , 1 d(F 2 ) A 2 dP 



(31) 



where 



d 2 ip | _ d 
dz 2 dw 



1 dip 
w dw 



Note that ( [28] ) and OTb are the equivalents of dAlb and JA31 ) in this 
geometry; see the Appendix for details. 

It is informative now to look in more detail at the torus equilib- 
ria found in the simulations. Once a conversion to a suitable cylin- 
drical coordinate system using the axis of symmetry of the equi- 
librium, quantities can be plotted in the (w, z) plane. In fig. [9] we 
see contours of the flux function ip as well as of F and P in the 
highest-helicity simulation of the set described in section [331 with 
Ai = 0.032; clearly the contours coincide, as in ( [301 ). We can also 
see that the cross-section of the flux tube is slightly elliptical close 
to its axis (the 'neutral line' where the poloidal field, or Bh, van- 
ishes) and a more complex shape further away; this is simply the 
tube's response to the torus geometry. 

In fig. [10] various quantities are plotted in the equatorial plane 
z = 0. Clearly, F(ip) is a smoother function than P(ip) \ it is possi- 
ble that fluctuations present during the formation of the torus have 
become 'frozen' into the equilibrium. Despite these fluctuations, 
the field is confirmed to satisfy ([311 to better than one part in 100 
throughout the volume. 



F 
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Figure 10. Plots of ip, F and P — P (multiplied for clarity by 5) against 
cylindrical radius vu along the line z = in the meridional plane, in the 
equilibrium illustrated in fig. [9] 



If the energies in the three components of the magnetic field 
(B w , Be/), Bz) are calculated as fractions of the total magnetic en- 
ergy, the values (0.172, 0.497, 0.331) are found. Whether the prox- 
imity of these fractions to 1/6, 1/2, 1/3 is coincidence is a matter 
of speculation at this stage. In the case of a straight flux tube, the 
fractional energies are 1/3, 1/3, 1/3 (see the Appendix). In addi- 
tion, it is found that if the quantity P ma g = B 2 /24tt is averaged 
over the entire volume, it is found to be within 2% of the average of 
the mean gas pressure difference P Q — P. That the average internal 
gas pressure is lower than the external pressure is a general result; 
the bubble cannot have a gas pressure excess confined by magnetic 
tension, an idea which is often found in the literature. It is however 
conceivable that a configuration exists where the gas pressure is 
higher than the external pressure at some region within the bubble, 
although its mean must still be lower than the external pressure. 



5 DISCUSSION 

The equilibria found in the simulations had low field strength, i.e. 
high plasma- f3. A low- f3 plasma would drive shocks into the exter- 
nal medium and could be numerically rather inconvenient. While 
there should be no fundamental difference between the structure of 
equilibria with f3 >> 1 and those with f3 w 1, since the pressure 
gradient force in both cases is able to balance the Lorentz force, 
an equilibrium with ^ < 1 should look rather different - it must 
be approximately force-free inside the bulk of the bubble; the pres- 
sure term drops out of ( [311 . In addition there must be a Lorentz 
force singularity on the boundary balanced by a discontinuity in 
gas pressure between the bubble interior and the external medium. 
Since the structure of such an equilibrium is constrained to a greater 
extent than that of a non-force-free equilibrium, it should be easier 
to construct such a field using analytic methods. Numerical investi- 
gation of this kind of equilibrium is left for the future. Note that the 
inhibition of the Rayleigh-Tayler instability at the bubble's leading 
edge might require a magnetic pressure comparable to the pressure 
of the ICM, meaning that we should expect disruption of a high- 
beta bubble. |PeYound ( [200i) discusses growth timescales and sta- 
bilising field strengths. 



5 In practice the 'discontinuity' should have some finite width, perhaps 
even non-negligible in comparison to the size of the bubble. 
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Also left for the future is a more thorough study of the he- 
licity and energy input from the AGN outflow. In section 12.31 we 
produced the greatest possible helicity by having the accretion disc 
and resultant jet threaded by a net flux, but it is not at all certain 
whether this should be the case in reality. The accretion disc can 
carry fl ux inwards from the sourc e cloud to the outflow launching 
region dSpruit & Uzdenskvll2005h but it is unclear what magnetic 
field geometry the source is likely to have. It is plausible that the 
net flux through the launching region fluctuates in time faster than 
a bubble inflates, so that the bubble contains regions with magnetic 
helicity of opposite signs. This will presumably mean that the first 
step once reconnection begins will be formation of localised equi- 
libria (as in the high-helicity case), followed by mutual annihilation 
and finally the formation of an extremely low-energy global equi- 
librium, if there is sufficient time. Another possibility is that the 
accretion does not have a significant net flux and the launching re- 
gion has mixed magnetic polarity. In this case, the helicity of the 
bubble is extremely low, even in small regions; it is not understood 
how the field should evolve in this case but it seems likely that re- 
connection should dissipate almost all of the magnetic energy on a 
rather short timescale. 

We have assumed here that there is no mixing of bubble mate- 
rial with the external medium. While the possibly of external mate- 
rial en tering the bubble by diffusive processes is ruled out by Pope 
j2010l) . it is possible that some external material becomes entrained 
in the bubble. First, there is the possibility of entrainment while the 
bubble is being inflated; the jet transfers momentum to the rest of 
the bubble in a potentially irregular m anner and external m aterial 
could easily become advected inwards. Crostonetal 1 d20Q8h for in- 
stance find evidence tha t in low-power radio sources (FR-I objects, 
iFanaroff & Rilevll 19741) . the pressure difference between the exter- 
nal medium and the radio-emitting component of the cavities is 
much more likely to be entrained external material than relativistic 
protons; the possibility of a strong magnetic field however remains. 
In addition, during the reconnection process there are blobs of ma- 
terial moving around at the Alfven speed (which may exceed the 
sound speed in the external medium) inside the bubble; it is plau- 
sible that some external matter becomes advected into the bubble. 
This material could become a non-magnetised 'island' inside the 
bubble or pass through a reconnection region and become magne- 
tised. Conversely, bubble material could leave the bubble as magne- 
tised islands, but would have to pass through a reconnection region 
in order to become magnetically disconnected from the rest of the 
bubble. 

It has also been assumed here that the gas in the bubble is at 
rest before reconnection begins. In reality, the gas is fired into the 
bubble at high speed and loses most of its kinetic energy in a shock 
but some will be left over, resulting in flows inside the bubble. Ne- 
glecting this motion simplified the analysis somewhat, but one can 
speculate on its possible effect. The fluid viscosity in such a dif- 
fuse medium is rather high so that such motions will be damped on 
a short timescale; to be more precise, r v isc = I 2 /v ~ / 2 /(Ac s ) 
where v, A and c s are the viscosity, mean free path and sound 
speed, respectively. With realistic parameters we find that the vis- 
cous timescale is comparable to the sound-crossing time. This is 
almost certainly less than the buoyant rise timescale calculated in 
dTTT) . However while these motions persist, they may prevent re- 
laxation to equilibrium if the kinetic energy is greater than the 
magnetic energy. This is unlikely in the case of a magnetically- 
accelerated outflow where only some fraction of the magnetic en- 
ergy is converted to kinetics 

As pointed out by Pfrom mer et aD <2005h . in the future it 



should be possible to use the Sunyaev-Ze rdovich (S-Z) effect to 
meas ure some of the u nknown parameter s dSunvaev & Zerdovichl 
1 19721 : for a review see Birkinshaw see lBasu et all 1201 d for 

the most recent results). Whereas the X-ray brightness is a line-of- 
sight integral, roughly speaking, of the square of the gas density, the 
S-Z intensity is a line-of-sight integral of the gas pressure. It may 
or may not be that we shall see dark cavities in the S-Z effect - it 
will be possible therefore to distinguish between a magnetically- 
dominated bubble ((3 <C 1) and a gas-pressure-dominated bubble 
(P ^> 1). In addition, by measuring the spectrum of the effect it is 
possible to distinguish between a non-relativistic and a relativistic 
plasma, which will resolve the issue of thermal gas vs. cosmic ray 
pressure. 



5.1 Other astrophysical contexts 

The equilibria found here show some similarity to those thought 
to reside in various kinds of non-convective star: upper-main- 
sequence, white dwarfs and neutron stars. A star which starts its 
life with some chaotic accretion/convection phase and a turbulent 
magnetic field undergoes the same kind of relaxation to equilib- 
rium once the convection dies away. The main difference between 
a star and a bubble is that in the former, gravity restricts motion in 
the radial direction. This provides extra stability, so that a greater 
range of stable equilibria is possible. For instance, it is possible in 
a star to have a much large r toroidal component than poloidal com- 
ponent (iBraithw aite 2009) whereas the two must be roughly equal 
in a bubble. Another way of thinking about this is that there are two 
forces - the pressure gradient — VP and gravity pg - which can be 
adjusted independently of each other to balance the Lorentz force, 
which also has two degrees of freedom (remember that V • B = 
removes one degree). 

There is also a striking similiarity between the twisted flux 
tubes found in the simulations and structures found in the iono- 
sphere of Venus. The Pioneer spacecraft measured the magnetic 
field vector as it orbited th e planet, frequently enco untering regions 
of strong magnetic field. [Russell & Elphid Jl979h found that the 
spacecraft was flying through twisted flux tubes. The tubes, about 
10km wide and with field strengths around 400 /xG, have a mag- 
netic pressure comparable to the external gas pressure, although 
it is not clear whether the gas pressure in the tubes is lower than 
or comparable to the magnetic pressure, i.e. whether the tubes have 
P <C 1 or p ~ 1. The surrounding ionosphere is much more weakly 
magnetised, with B « 20 ^ G. Similar tubes have since been found 
in the ionospheres of Mars dCloutier et al.ll 1999b and Titan. From 
the considerations in section H~T1 it should not be surprising that 
the tubes are twisted, as that state is in some sense the 'natural' 
condition of a flux tube. However, any flux tube needs to be held 
at the ends, either by a torque (in the case of a twisted tube with 
2^ax = ^h) or by a tension (in the case of an untwisted tube). 



6 CONCLUSIONS 

I have considered the evolution of the magnetic field inside AGN- 
inflated bubbles which are observed as dark cavities in X-ray im- 
ages of galaxy clusters. It is found that the magnetic field under- 
goes relaxation to a global- scale equilibrium filling the entire bub- 
ble, consisting of twisted flux tube(s) arranged in some pattern. The 
relaxation process inevitably involves magnetic reconnection - the 
reconnection regions could provide energetic synchrotron-emitting 
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particles via X-point and Fermi acceleration (see e.g. Parker 19571 : 
lMilleretal.ll 19971) . 

The timescale on which this relaxation takes place, or in other 
words the stage during this relaxation we are likely to observe, de- 
pends crucially on various parameters: the magnetic field strength, 
mass density, Lorentz factor and size of the outflow as well as the 
properties of the ambient intra-cluster medium into which the bub- 
ble expands. Given the uncertainly in these parameters, it is impos- 
sible at this stage to distinguish between the following eventualities 
(see section [2T2t . In the following, the radius of the bubble and the 
dominant length scale of its magnetic field structure are r and / 
respectively; the Alfven speed is va = B/ ^Anp and there is a 
reconnection timescale T rec = 1/{olva) where a « 0.1 is the re- 
connection speed parameter. 

(i) The AGN outflow is weakly magnetised and little reconnec- 
tion occurs; the observed field is small-scale and evolves passively 
in response to the bubble's interaction with the ICM. Measurement 
of the relevant parameters would show that / <C r and T rec > T age 
where T age is the age of the bubble. 

(ii) The AGN outflow is strongly magnetised and the magnetic 
field relaxes towards a global equilibrium. However, because the 
helicity of the field is low, the bulk of the magnetic energy is dis- 
sipated and no global equilibrium is reached. At the time of obser- 
vation the reconnection is still ongoing and T rec ~ T age . The field 
may consist of local-equilibrium flux tubes of size I < r. 

(iii) The AGN outflow is strongly magnetised and has high he- 
licity so that reconnection proceeds quickly and a global equilib- 
rium is reached: / « r and T rec > T age . The field may consist 
of large-scale twisted flux tube(s) arranged in figure-of-eight pat- 
terns or as a single torus configuration, similar to the spheromak 
shape found in laboratory experiments. In this case, the magnetic 
field will give the bubble some rigidity, helping keep it intact as it 
moves through the ICM. 

To illustrate this with plausible parameters, if we measure a den- 
sity 10 -5 ra p g cm -3 and field strength 20/iG in a bubble of radius 
lOkpc then T rec ~ (l/r) 7Myr; if the bubble is older than 7Myr 
then we have the global-equilibrium case (iii). To reach this situa- 
tion the AGN outflow must have had high helicity; this is likely if 
the accretion disc is fed material with a consistent net flux. Fluctu- 
ating or vanishing net flux through the accretion disc will result in 
case (ii) even if the outflow is strongly magnetised. 

During reconnection to equilibrium, the shape of the bubble 
may change in response to plasma flow inside the bubble on the 
order of the Alfven speed. However, if the density of the bubble is 
much less than the density of the surrounding ICM, the effect on 
the shape of the bubble will be rather modest. 

Finally, it is shown that the difference in gas pressure between 
a bubble and its surroundings is equal to one third of the magnetic 
energy density, i.e. the magnetic field produces an 'isotropic mag- 
netic pressure' P mag = B 2 /24tt — P — P{ where the subscripts o 
and i denote pressure outside and inside the magnetised volume. In 
this and other contexts this is more useful than the B 2 /8tt which 
is more common in the literature; this is a general feature of three- 
dimensional problems. 
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APPENDIX A: 
TUBE 



STRUCTURE AND TENSION OF A FLUX 



Consider a straight flux tube, uniform along its length, with some 
arbitrary cross-section. The field component parallel to the direc- 
tion of the tube is P ax , whilst Bh represents the other two com- 
ponents which can be thought of as field lines in the plane of the 
cross-section0 We can see from the argument in section 01 that in 
equilibrium the pressure P is constant on lines of Bh - Furthermore, 
the axial component of the Lorentz force must vanish in a tube 
which is uniform along its length, since there is no axial pressure 
gradient to balance it. From this it can be shown that: 



B h -VP ax = 0, 



(Al) 



which in other words means that P ax must be constant along lines 
of Bh- We can now describe the magnetic and pressure fields in 
terms of the axial component A ax of the vector potential: 

B h = V x (A ax e ax ); P = P(A ax ); P ax = P ax (A ax ); (A2) 

where e ax is the axial unit vector. Equating the Lorentz force to the 
pressure gradient gives 

dP 



V 2 A 



+ 



1 d(Bi 



+ 4:71- 



0. 



(A3) 



2 dA ax dA a 

The three equations dAll ) to ([A3} are the equivalent of equations 
(f28t to (THl in this geometry. 

In addition to this, there is a good reason to believe that 
tubes will have a circular cross-section. Using cylindrical coordi- 
nates (r, 0, z), imagine perturbing a tube with circular cross sec- 
tion, radius ro and field Bh = e^P^ofV) with a displacement 
field £ = e r rx(0) where r, e r , cj) and are the radial and az- 
imuthal coordinates and unit vectors and x is an arbitrary function. 
If r = r + £ = (1 + x)r then dr' = (1 + x)dr and the magnetic 
energy per unit length of the Ph component is 

2tt 



hr r 

O7r «/(/) = Jr = 



r'dr'dcj) [B[ 



, + B 2 r ] 



(A4) 



r(j>=1ir p r — r 

8vr L = J r = 



(1 + x) 2 r drdcj) 



Bs 



1 + x y 



+ B 2 



where B^ = B^o/fl + x) follows from flux freezing. Clearly 
the energy in the <j) component of the field is unchanged by the 
perturbation but for any function x which is not independent of 0, 
i.e. that gives the tube a non-circular cross-section, there is new 
energy in the B r component which was absent before, meaning 
the total energy has increased. Circular tubes therefore represent an 
energy minimum. The first term in the equilibrium condition dA3l) 
simplifies to (l/r)(d/dr)(r dA ax /dr). 

Now consider such a tube with circular cross-section of ra- 
dius a, length / and volume V = no 2 1 containing r.m.s. axial and 
azimuthal field components P ax and Bh . First I examine the ener- 
getics and stability of the tube to various pertubations; later I look 
at the boundaries at either end of the tube. 

The axial and hoop fluxes $ ax and <E>h of the magnetic field 
are given by 



2 d 
7VCL B a 



and 



k h <$>h = alB h 



(A5) 



where /c ax and kh are dimensionless factors of order unity which 
correct for the different averaging required when calculating energy 
and fluxQ The fluxes are conserved on a dynamical timescale, as 
is the magnetic helicit)[f| of the tube which can be expressed as 
H ~ $ ax $ h - The magnetic energy of the tube is given by 



E 



V_ 

-(- 

8 \tvV 
1 f V 



(5 ax + B h ) 
2 k 2 $ 2 



+ \k& 2 h 



:k 2 <$> 2 



KCL l2l 2 



(A6) 
(A7) 

(A8) 



Furthermore, as the two fluxes are constant during dynamic adjust- 
ments the derivatives w.r.t. V at constant / and a are 



dE_ 

dv 

8E 
dV 



]2i2 ^2 

8tvV 2 

k 2 $ 2 
87r 3 a 4 ~ 



8tt 



8V 2 



^ax 
87T 



8tt 



(A9) 



(A10) 



respectively, using dA5b and specifying that the derivatives repre- 
sent homogeneous expansion where k ax and kh are constant. Now 
defining a magnetic pressure given by P ma g = —dE/dV, we see 
from JA91 ) that a tube of fixed length provides a positive pressure in 
the lateral direction. Therefore once a dynamical equilibrium has 
been reached the average thermal pressure in the tube Pi should be 
lower than that in the surroundings P Q by a quantity 



Pi 



Pax 

8tt 



(All) 



The same must be true in the axial direction, as the thermal pressure 
acts the same in both directions and so therefore must the Lorentz 
force. This means that the derivatives ( IA91 ) and dAlOl ) are equal, 
giving 



Hh — ^ax 



(A12) 



6 The suffices ax and h stand for 'axial' and 'hoop'. The reader will see 
that finding a consistent terminology for both this case and the case of the 
axisymmetric equilibria described in section l4~2l is less than straightforward. 
In the literature, normally z and 4> are used for the axial and hoop directions 
respectively. However, when a tube is connected into a circular loop the 
axial direction becomes 'azimuthal' and is denoted by <f>, and z is now the 
axis of the loop. 



7 Instead of using averages one could consider a particular flux tube, for 
instance with functions P ax (r) = const and £>h( r ) — Bor/a, but this 
provides no extra insight. 

8 Helicity is gauge independent if the domain is bounded by magnetic sur- 
faces and/or periodic conditions. Here, we have the former on the sides of 
the tube and can assume the latter for the ends. 
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It is also possible to imagine a tube which has reached equilibrium 
in the lateral direction so that relation dAllb is satisfied, but which 
for reasons to do with whatever it is anchored to at the ends, is not 
in equilibrium in the axial direction. In other words, the tube has a 
net tension or pressure along its length. This tension is calculated 
thus: 

Tva 2 = (2Bl x - B 2 )a 2 /8, (A13) 

using dAlOb and dAllb . The term P — Pi comes from the fact that 
in stretching the tube, work P Q dV must be done against the exter- 
nal medium while the internal gas does work P\dV . Alternatively, 
to avoid doing PdV work the tube may be stretched at constant 
volume, in which case the tension T = (^f)v, which gives the 
same result. Of course, the result cannot depend on the change in 
a during the stretching, because lateral force balance has already 
been assumed and changes in a are energetically neutral. Also note 
that the equilibrium condition (IA121 ) corresponds to vanishing ten- 
sion T0 Finally, note that as assumed in © we have 

P mag = (A14) 
where B is the total magnetic field given by B 2 = B 2 ^ + B\. 

This paper has been typeset from a TgX/ ETjnX file prepared by the 
author. 



dV 1 + ° 



Pi 
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9 IWeissI dl964h and authors of various later works erroneously state the 
tension of a flux tube as (£>ax — B^)a 2 /8, neglecting the PdV work. 
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